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Abstract 

Background:Tw'\sU and Twist2 are highly homologous bHLH transcription factors that exhibit extensive highly overlapping 
expression profiles during development. While both proteins have been shown to inhibit osteogenesis, only Twistl 
haploinsufficiency is associated with the premature synostosis of cranial sutures in mice and humans. On the other hand, 
biallelic Twist2 deficiency causes only a focal facial dermal dysplasia syndrome or additional cachexia and perinatal lethality 
in certain mouse strains. It is unclear how these proteins cooperate to synergistically regulate bone formation. 

Methods: Twistl floxed mice {Twist 7 f/f ) were bred with Twist2-Cre knock-in mice {Twist2 Cre/+ ) to generate Twistl and Twist2 
haploinsufficient mice {Twist 7 f/+ ; Twist2 c ' e/+ ). X-radiography, micro-CT scans, alcian blue/alizarin red staining, trap staining, 
BrdU labeling, immunohistochemistry, in situ hybridizations, real-time PCR and dual luciferase assay were employed to 
investigate the overall skeletal defects and the bone-associated molecular and cellular changes of Tw/sf 7 f/+ ,Tw/sf2 Cre/+ mice. 

Results: Tw\sf\ and Twist2 haploinsufficient mice did not present with premature ossification and craniosynostosis; instead 
they displayed reduced bone formation, impaired proliferation and differentiation of osteoprogenitors. These mice 
exhibited decreased expressions of Fgf2 and Fgfrl-4 in bone, resulting in a down-regulation of FGF signaling. Furthermore, 
in vitro studies indicated that both Twistl and Twist2 stimulated 4.9 kb Fgfr2 promoter activity in the presence of El 2, a 
Twist binding partner. 

Conclusion: These data demonstrated that Twistl- and 7VWsf2-haploinsufficiency caused reduced bone formation due to 
compromised FGF signaling. 
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Introduction 

Mammalian Twistl and Twist2 are two members of the Twist 
subfamily of the basic-helix-loop-helix (bHLH) transcription 
factors that have been highly conserved during evolution [1]. In 
Drosophila, a single Twist gene, DTwist, is essential for embryonic 
gastrulation and mesodermal formation [2,3]. Mouse Twistl was 
identified by its high homology with DTwist [4,5], while Twist2, 
originally called "Dermol", was discovered by a yeast-two-hybrid 
screen using the ubiquitous bHLH protein El 2 as bait [6]. The 
expression patterns of Twistl and Twist2 show an extensive 
overlap during mouse embryonic development [6], and their 
encoded proteins exhibit a high degree (up to 98%) of sequence 
similarity [7] . Both proteins perform various functions by forming 
either homodimers or heterodimers with bHLH E proteins (El 2/ 
E47) that bind to DNA canonical regulatory sequences called "E- 
boxes" (CANNTG) in the promoter region of target genes [8]. 

In humans, mutations in the T1VIST1 gene are associated with 
Saethre-Chotzen Syndrome (SCS), which is an autosomal 



dominant disorder characterized by craniosynostosis, brachydac- 
tyly, soft tissue syndactyly and facial dysmorphism [9] . The skeletal 
phenotype of Twist 1 -heterozygous mouse consistently resembles 
that of human SCS with premature fusion of the cranial suture 
[9,10]. As mouse embryonic development progresses, the Twistl 
expression declines in the developing bones of the skull [1 1] . In 
addition, Twistl overexpression was found to inhibit osteoblast 
differentiation in vitro and in vivo [12,13,14]. Together, these 
observations suggest that Twistl negatively regulates osteoblast 
differentiation and bone formation. 

Various molecular mechanisms may be responsible for the 
inhibitory role of Twistl in osteoblast differentiation. Twistl may 
modulate FGF signaling, especially Fgfr2 expression in cranial 
suture development [15,16,17] or it may directly bind to and 
inhibit the transactivation function of Runx2, a master regulator of 
osteogenesis [1 1], In addition, Twistl might indirectly regulate the 
Runx2 expression through modulating FGFR2 expression as shown 
in the ex vivo cultured primary osteoblasts isolated from human 
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Table 1. Primers used for real-time PCR. 





Gene 


Forward primer (5 -3') 


Reverse primer (5 -3') 


Fgf2 


agcggctctactgcaagaac 


gccgtccatcttccttcata 


Fgfrt 


ttcagtggctgaagcacatc 


gcagagtgatgggagagtcc 


Fgfr2 


gtgcttggcgggtaattcta 


gatgactgtcaccaccatgc 


Fgfrt 


ccaccttcaagcagttggtag 


gggtgaacaccgagtcatct 


Fgfr4 


aagtcatccgtggccactac 


cgagaggcaggtctagattca 


Erm 


gatgatgcctgaaagccagt 


gggaccccatgttcatacag 


Pea3 


cctatgactcccccagacaa 


cctccctgaggagatgtgaa 
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SCS patients [18]. Finally, it is possible that Twistl inhibits 
osteoblast apoptosis via the suppression of TNF-oc expression [19]. 

Twist2 has been shown to have an inhibitory function similar to 
that of Twistl in bone formation [11]. While recessive TWIST2 
mutations in humans and its inactivation in mice result in a focal 
facial dermal dysplasia (FFDD) syndrome, there is no Twist2- 
deficient skeletal phenotype [20]. The phenotypic difference 
between the Twistl- and Twist2-deficient subjects is indeed 
intriguing when viewed in the context of their significandy 
overlapping expression patterns in vivo [6] and their similar 
functions in bone formation [1 1]. Thus, it is largely unknown how 
Twistl and Twist2 synergistically regulate bone formation and 
what molecular mechanism is involved. 

In this study, we generated a compound Twistl- and Twist2- 
haploinsufficient animal model, Twistr 0 " + ; Twist2 Cre/+ mice, by 
crossing Twistl floxed mice with Twist2 Cre knock-in mice. Thus, 
the compound Twistr" + ; Twist2 c " !/+ mice had one allele of 
Twist2 replaced by the cre recombinase and one allele of Twistl 
deleted specifically in tissues where the Twist2 gene was expressed. 
To our surprise, we found that the TwHstF"" + ; Twist2 c " /+ mice 
exhibited reduced bone formation and impaired proliferation and 
differentiation of osteoprogenitors. The skeletal abnormalities were 
associated with reduced FGF signaling as a consequence of the 
decreased expressions of Fgf2 and Fgfrs. 

Materials and Methods 

Ethical Approval 

All animal protocols were approved by the Institutional Animal 
Care and Use Committee (IACUC) at Texas A&M University 
Baylor College of Dentistry. IACUC has specifically given ethical 
approval for all the procedures in this study. 

Generation of Twistl flox/+ ; Twist2 Cre/+ 

Twistl floxed mice (Twist P 1 ™^*) were maintained on a C57/ 
1 29 mixed genetic background [2 1] and Twist2 Cre Knock-in mice 
(Twistr"'-*) were on a 129 genetic background [22]. The 
Twist P 1 "*^ 0 * mice were bred with Twist2" e/+ mice to generate 
compound TwistF 7 ° x/+ ; Twist2 cre/+ mice. The Twist2 c " /+ mice 
carry a Cre recombinase that replaces one allele of the Twist2 gene 
[22]; therefore, the floxed Twistl allele is deleted in the tissues 
where the Twist2 gene is active. The Twistl floxed mice and 
Twist2 c " /+ mice were genotyped as described previously [21,22]. 
In this study, we analyzed the skeletal phenotype of 6-8 day-old 
TwistP 1 ""^; Twist2 Cre/+ mice and used age-matched Twistr* 1 ™ 1 *; 
Twist2 +/+ as control mice (counting the day of birth as day 0). 



Alcian blue/alizarin red staining of the skeleton 

Alcian blue/ alizarin red staining was performed to analyze the 
overall skeletal and mineralization defects, as described previously 
[23]. Briefly, 6-day-old Twist Twist2 cre/+ mice and control 
mice were sacrificed, skinned, eviscerated and fixed for three days 
in 95% ethanol. They were then stained with alcian blue for 
cartilage and alizarin red for bone visualization. 

Plain X-ray radiography and high-resolution 
microcomputed tomography (u-CT) 

The femurs and tibiae from 6-day-old TwistV 7 ""^; Twist2 cre/+ 
mice and control mice were dissected free of the skeletal muscles 
and fixed in 70% ethanol. For plain X-ray radiography, the 
femurs were analyzed with a Faxitron MX-20 specimen radiog- 
raphy system (Faxitron X-ray Corp., Buffalo Grove, IL) as 
described previously [24] . For the high-resolution u-CT analyses, 
the tibiae were scanned at 3.5-um resolution using a U.-CT35 
imaging system (Scanco Medical, Basserdorf, Switzerland), as 
previously described [24]. The trabecular bone was analyzed at a 
threshold of 160 in 20 sections underneath the growth plate. 

Histology, immunohistochemistry and in situ 
hybridization 

For histologic analysis, the bone specimens from 6-day-old mice 
were fixed in freshly prepared 4% paraformaldehyde, decalcified 
in 10% EDTA with 0.5% paraformaldehyde, and embedded in 
paraffin using standard procedures [25] . Serial 7-um sections were 
cut and mounted on silane-coated slides. The sections were then 
used for Hematoxylin and Eosin (H&E) staining, Tatrate-resistant 
alkaline phosphatase (TRAP) staining, immunohistochemistry or 
in situ hybridization, as described previously [24,25] . 

The following antibodies were used for the immunohistochem- 
ical analyses: anti-Osterix (Santa Cruz Biotechnology, Dallas, TX, 
USA; 1:400), rabbit anti-biglycan antibody "LF-159" (a gift from 
Dr. Larry Fisher, National Institutes of Health, Bethesda, MD; 
1:1000), anti-Fgfr2 polyclonal antibody (Abeam, Cambridge, MA, 
USA; 1:5000), and anti-phospho-p44/42 MAPK (Erkl/2) anti- 
body (Cell Signaling Technology, Danvers, MA, USA; 1:100). All 
the immunohistochemical experiments were performed with a 3, 
3-diaminobenzidine kit (Vector Laboratories, Burlingame, CA) 
according to the manufacturer's instructions. 

For the in situ hybridization, the RNA probes for dentin matrix 
protein (Dmpl), osteocalcin (Ocm) and alkaline phosphatase (Alp) 
were labeled with digoxigenin (DIG) using a RNA Labeling Kit 
(Roche, Indianapolis, IN, USA) according to the manufacturer's 
instructions. DIG-labeled RNA probes were detected by an 
enzyme-linked immunoassay with a specific anti-DIG-alkaline 
phosphatase antibody conjugate (Roche, Indianapolis, IN) and a 
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Twist2 Cre/+ Twist2 Cre/+ Twist2 Cre/+ 

Figure 1. Reduced bone formation in Twist f ox/+ ; Twist2 Cre/+ mice. (A) Skeletons of 6-day-old control (left) and Twist f"^*; Twist2 Cre/+ (right) 
mice stained with alcian blue (cartilage) and alizarin red (bone). The skeleton of the Twist1 ,lox/+ ; Twist2 Cre/+ mouse is remarkably smaller. (B) Alcian 
blue- and alizarin red-stained skull from 6-day-old Twist1 flox/+ ; Twist2 Cre/+ mice (right) showed delayed fusion of interfrontal suture and open posterior 
fontanel (arrows), compared with the control mice (left). (C) Alcian blue- and alizarin red-stained hind foot of 6-day-old control (left) and Twist1 flox/+ ; 
Twist2 Cre/+ (right) mice. Note the delayed ossification in metatarsals (mt) and phalanges (pi), and an additional toe (arrow) originating from the same 
(or duplicated) metatarsal as the hallux in Twistf° x/+ ; Twist2 Cre/+ mice. (D) Plain X-radiography of the tibiae from 6-day-old control (left) and Twistf ox/+ ; 
Twist2 Cre/+ mice (right). The Twistf ox/+ ; Twist2 Cre/+ mice had shorter tibiae and reduced radiopacity, compared to the control mice. (E) 
Representative three-dimensional u-CT images of tibiae from 6-day-old control (left) and Twist1 flox/+ ; Twist2 Cre/+ (right) mice. The Twist1 flox/+ ; 
Twist2 Cre/+ mice showed reduced trabecular (arrowheads) and cortical bones (arrows). (F-H) Quantitative |x-CT data showing that the 6-day-old 
Twist f ox/+ ; Twist2 Cre/+ mice had a significant decrease in the ratio of bone volume (BV)/total volume (TV) (F) and in apparent bone density 
(G), compared to the control mice (n = 6, P<0.001). The Twist1 ,lox/+ ; Twist2 Cre/+ mice also presented reduced material density although no 
statistically significant difference was observed (H). 
doi:10.1371/journal.pone.0099331.g001 
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Control Twistlfl ox/+ ; Twist2 CreA Control Twist l^ ox/+ ; Twist2 Cre/+ 




Figure 2. Histological examination of Twist f° x/+ ; Twist2 Cre/+ mice. (A-B) Femur sections of 6-day-old control and Twistf ox/+ ; Twist2 Cre/+ mice 
were stained with H&E. The Twist1 flox/+ ; Twist2 Cre/+ mice displayed reduced metaphyseal trabecular bone (A, red arrows) and a decreased thickness of 
the periosteum (B, blue arrows) and cortical bone (B, red arrows). (C) TRAP staining of femur sections of 6-day-old control and Twist1 flox/+ ; Twist2 Cre/+ 
mice. Note that the osteoclasts (red arrows) appeared to be similar in size and distribution in the control and Twist1 flox/+ ; Twist2 Cre/+ mice. The 
osteoclast densities were 0.55±0.06/0.01 mm 2 in the controls (n = 5) and 0.60 ±0.02/0.01 mm 2 in the Twistf ox/ +; Twist2 Cre/+ mice (n = 5, P>0.05). (D-F) 
In situ hybridization analyses (signal in blue) of the transcripts of Alp (D), Ocn (E) and Dmpl (F) in the femurs of one-week-old control and Twist1 flox/+ ; 
Twist2 Cre/+ mice. (G, H) Immunohistochemical analyses (signal in brown) of the osterix (G) and biglycan (H) protein levels in the femurs of the 6-day- 
old control and Twistl" ox/+ ; Twist2 Cre/+ mice. Scale bar= 100 \im. 
doi:1 0.1 371 /journal.pone.0099331 .g002 



VECTOR NBT/BCIP alkaline phosphatase substrate (Vector 
Laboratories, Burlingame, CA), which produced a blue color 
indicating positive signals. Methyl green was used for counter- 
staining. 

Cell proliferation assay 

Six-day-old TwistF 7ax/+ ; Twist2 c " /+ mice and control mice were 
intraperitoneally injected with 5-bromo-2'-deoxyuridine (BrdU) 
( 1 0 mg/ 1 00 g body weight) at 24 hours and then two hours before 
sacrifice. The long bones were collected and processed for paraffin 
sectioning as described above. The incorporated BrdU was 
detected with a BrdU staining kit (Invitrogen, Camarillo, CA, 
USA) according to the manufacturer's instructions. The BrdU- 
positive nuclei and total nuclei were counted in the metaphysis 
region (trabecular bone) as well as in the midshaft (cortical bone). 
The data represented the percentage of BrdU-positive nuclei from 
four individual animals each group. 

Quantitative Real-time PCR 

Total RNA was extracted using Trizol (Invitrogen, Camarillo, 
C A, USA) from the long bones of 8-day-old Twist 1 / 2 dHet mice 



and control mice and reverse-transcribed into cDNA with a 
Reverse Transcription Kit (QIAGEN, Germantown, MD, USA). 
Quantitative real-time PCR (qPCR) was performed using a Go 
Tag qPCR Master Mix System (Promega, Madison, WI, USA). 
Mouse 18S rRNA primers (PPM57735E-200, QIAGEN, German- 
town, MD, USA) were used for normalization. The primers for 
Alp, Ocn, bone sialoprotein (Bsp), osterix (Osx), Dmpl and Runx2 
were reported elsewhere [26,27]. The primers for Fgf2, Fgfrl, 
Fgfr2, Fgfr3, Fgfr4, Erm (EtvS) and Pea.3 [Etv4) are listed in Table 1. 
The primers for Fgfr2 were specifically designed to amplify the 
mesenchymal isoform of Fgfr2 (Fgfi2-IIIc). All experiments were 
performed in triplicate on three animals. The data were analyzed 
using the 2 AAGT method as described previously [28] . 

Cell culture, constructs and promoter luciferase assay 

The C3H10T1/2 mesenchymal cells and MC3T3-E1 preosteo- 
blast cells were cultured as described previously [29]. A 4.9 kb- 
Fgfr2 promoter luciferase construct and expression constructs for 
Twist 1 and El 2 have been reported elsewhere [30]. An expression 
construct for Twist2 was generated by cloning Twist2 cDNA into 
the BamHI and EcoRI sites of pcDNA3 vector (Invitrogen). The 
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Figure 3. Quantitative real-time PCR analyses of osteoblast differentiation markers. Real-time PCR was performed with total RNA isolated 
from the long bones of the one-week-old control and Twist1 flox/+ ; Twist2 Cre/+ mice. The expressions of key transcription factors associated with 
osteoblast differentiation (Runx2 and osterix), osteoblast markers (Alp, Ocn and Bsp) and osteocyte marker [Dmp1) were reduced in Twistl" 0 "'*; 
Twist2 Cre/+ mice. The mRNA levels in the control mice were set as one, and the mRNA levels of Twist1 flox/+ ; Twist2 Cre/+ mice were expressed as folds of 
those in the control mice. The data represented three analyses (n = 3) for each group. 
doi:1 0.1 371 /journal.pone.0099331 .g003 



promoter luciferase assay was performed as described previously 
[24,30]. Briefly, C3H10T1/2 cells and MC3T3-E1 cells were 
plated in 24-well plates at a density of 3 x 10 4 cells per well. Then 
the cells were transiently transfected with 0.1 |ig of 4.9 kb-Fgfr2 
promoter luciferase construct, together with 0.4 u.g of various 
constructs expressing Twist 1, Twist2 or E12 using FuGENE 6 
Transfection Reagent (Roche, Indianapolis, IN, USA). The total 
amounts of transfected DNA were balanced by the addition of an 
empty vector (pcDNA3). A Renilla luciferase expression construct 
was co-transfected as an internal control to monitor the 
transfection efficiency. Forty-eight hours later, the transfected 
cells were analyzed using a dual luciferase reporter assay system 
(Promega, Madison, WI, USA). The luciferase activities were 
normalized by the control. All experiments were carried out in 
triplicate and repeated three times. 

Statistical analysis 

The statistical analyses were performed with a one-way 
ANOVA for a multiple group comparison and the Student's t- 
test for a two-group comparison. If significant differences were 
found with the one-way ANOVA, the Student's /-test was used to 
determine which groups were significandy different from the 
others. The quantified results were expressed as mean ± standard 
deviation (SD). P<0.05 was considered to be statistically 
significant. 

Results 

Generation of Twistl- and Twist2-haploinsufficient mice 

Twistl^" x/+ ; Twist2 Cre/+ mice were generated by breeding Twistl 
floxed mice (Twist ]fi" xl fl'"^ with Twist2 Cre Knock- in mice 
(Twist2 c " /+ ), in which one allele of Twist2 was replaced by the 
Cre recombinase. Twist2-Cre is active in the condensed mesen- 
chyme that will later produce the chondrocytes and osteoblasts 
[22]; therefore, in the compound Twist f" x/+ ; Twist2 Cre/+ mice, the 
floxed Twistl allele is deleted in the osteoblasts and their 
precursors. Consequently, both Twistl and Twist2 were haploin- 
sufficient in the osteoblast lineage of the Twistl^" x/+ ; Twist2 Cre/+ 
mice. Real-time PCR confirmed that the mRNA levels of both 
Twistl and Twist2 were about three folds less in the Twist f" x/+ ; 
Twist2 Cre/+ mice than in the control mice (Figure SI). Although the 
floxed Twistl allele was conditionally deleted, the compound 
Twist f 1 "^; Twist2 Cre/+ mice had a phenotype similar to that of the 
mice completely heterozygous for both the Twistl and Twist2 



genes [31]. Most of them died within two weeks after birth but a 
few of them survived to adulthood and were fertile. 

Reduced bone formation in Twistl flox/+ ; Twist2 Cre/+ mice 

We first examined the overall skeletal structures of Twistl^ ox/+ ; 
Twist2 " + mice in Figure 1. The alcian blue/alizarin red staining 
showed that TwistV 1 "'*; Twist2 c " /+ mice had a much smaller 
skeleton with delayed fusion of the interfrontal suture, open 
posterior fontanelles and delayed ossification in the metatarsals 
and phalanges (Figs. 1A-C). In addition, the Twutf"' /+ ; Twist2 Cre/+ 
mice developed an extra toe close to the hallux (Fig. 1C), the 
hallmark of Twistl heterozygous mice [32]. Plain X-radiography 
showed that the Twistf" xU ; Twist2 c " /+ mice had reduced 
radiopacity in the tibiae compared to the control mice (Fig. ID). 
Micro-CT images further confirmed that the Twistl^° x/+ ; 
Twist2 c " /+ mice had reduced trabecular bone and decreased 
cortical bone thickness (Fig. IE). The quantitative analyses 
revealed a significant decrease in trabecular bone volume versus 
total volume (BV/TV) and in apparent bone density in the 
TwutP l " x/+ ; Twist2 Cn/+ mice (Figs. 1F-G). The material density 
was also slightly reduced in the Twistl^ hx/+ ; Twist2 CTe/+ mice 
although the difference in this parameter was not statistically 
significant compared to the control mice (Fig. 1H). These data 
demonstrated that bone formation was inhibited in the 
Twistf' ,x/+ ; Twist2 Cre/+ mice. 

Defects in osteoblast differentiation in Twistl flox/+ ; 
Twist2 Cre/+ mice 

We carried out a series of histological and molecular analyses to 
determine whether there were any abnormalities in osteoclast 
numbers and/or osteoblast differentiation. First, H&E staining 
confirmed that the Twutl-^" x/+ ; Twist2 c " /+ mice formed less 
trabecular bone and thinner periosteum and cortical bone 
(Figs. 2A-B). Second, TRAP staining showed that the distribution 
and size of the osteoclasts in the Twist f 7 '"^; Twist2 c " l/+ mice were 
similar to those of the control mice (Fig. 2C). Although the 
osteoclast density was slightly increased in the TwistP 10 "^; 
Twist2 c " /+ mice, no significant difference was observed between 
two groups. Third, we examined the expression levels of the 
osteoblast differentiation markers by in situ hybridization, immu- 
nohistochemistry and real-time PCR. These methods revealed that 
the levels of the osteoblast differentiation markers Alp (Figs. 2D 
and 3), Ocn (Figs. 2E and 3), biglycan (Fig. 2H) and Bsp (Fig. 3) 
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Figure 4. Cell proliferation was reduced in the Twist f ox/+ ; Twist2 Cre/+ mice. (A-B) BrdU immunohistochemical staining of the femur sections 
of 7-day-old control and Twist1 fl ° x/+ ; Twist2 Cre/+ mice. The BrdU-positive cells (signal in brown) were counted in a 100-um zone of the metaphysis, 
demarcated by the chondro-osseous junction and the marked line (A), and in the femoral diaphysis (B). The osteoblast/osteoprogenitor proliferation 
was reduced in both the metaphysis (C) and diaphysis (D) in the Twist 1 flox/+ ; Twist2 Cre/+ mice, compared to the control mice (n = 4, *P<0.05). The data 
were expressed as the percentage of BrdU-positive nuclei versus total nuclei. Scale bars: 100 urn. 
doi:1 0.1 371 /journal.pone.0099331 .g004 



were sharply reduced in the Twistl^°" /+ ; Twist2 Cre/+ mice 
compared to the control mice. In addition, the osteocyte marker 
Dmpl was also dramatically decreased (Fig. 2F and 3). Further- 
more, the expression levels of Runx2 and Osx, two key 
transcription factors essential for osteoblast differentiation and 
bone formation, were remarkably reduced in the TwistP 1 "*'*; 
Twist2 Cre/+ mice, compared to the control mice (Figs. 2G and 3). 
Taken together, these findings supported the hypothesis that the 
reduced bone formation resulted from defects in osteoblast 
differentiation rather than abnormal osteoclast activities. 



Reduced cell proliferation in Twist f° x/+ ; Twist2 Cre/+ mice 

Since the Twistl^° x/+ ; Twist2 Cre/+ mice had a reduced periosteal 
layer compared to the control mice (Fig. 2B), we performed BrdU 
labeling to determine whether the proliferation of osteoblasts and 
their progenitors was affected. We noticed that the BrdU-positive 
cells in the area of the metaphysis (Figs. 4A), as well as in the mid 
diaphyseal periosteum and cortical bone (Fig. 4D), were signifi- 
cantly reduced in the Twistr x/+ \ Twist2 c " /+ mice. However, no 
difference in the osteoblast apoptosis was observed (Figure S2). 
These findings suggested that the reduced cell proliferation might 
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Figure 5. Reduced FGF signaling in Twistl flox/+ ; Twist2 Cre/+ mice. (A) Quantitative real-time PCR was performed to analyze the mRNA levels of 
Fgf2, Fgfrl, Fgfr2, Fgfr3 and Fgfr4, as well as two effector molecules, Erm and Pea3, using total RNA isolated from the long bones of 8-day-old control 
and Twistf ox/+ ; Twist2 Cre/+ mice. The mRNA levels in the control mice were set as one, and the mRNA levels of Twistl"° x/+ ; Twist2 Cre/+ mice were 
expressed as folds of those in the control mice. The data represented three analyses (n = 3) for each group. (B) Immunohistochemistry showed that 
the level of the Fgfr2 protein (signal in brown) was reduced in the femurs of the Twist1 flox/+ ; Twist2 Cre/+ mice (right), compared to the control mice 
(left). (C) Immunohistochemistry showed that there was less phospho-Erk1/2 (signal in brown) than in the femurs of Twist1 flox/+ ; Twist2 Cre/ * mice, 
compared to the control mice. Scale bars: 100 um. 
doi:1 0.1 371 /journal.pone.0099331 .g005 



also contribute to the reduced bone formation in the TwistF 7ox/+ ; 
Twut2 c " /+ mice. 

Reduced FGF signaling in Twist1 flox/+ ; Twist2 Cre/+ mice 

A variety of studies have indicated interactions between Twist 1 
and FGF signaling although the outcome of such interactions 
appears to be context dependent [15,17,18,33]. We analyzed the 
components of the FGF signaling pathway in the Twistf^° xl+ ; 
Twut2 C " /+ mice. Quantitative real-time PCR revealed a substan- 
tial decrease in the transcript levels of Fgf2 and Fgfrsl, 2, 3 and 4 in 



the Twist f° x/+ ; Twist2 c " /+ mice (Fig. 5A). Immunohistochemistry 
further confirmed that the Fgfr2 protein was markedly reduced in 
the long bones of the TwistP 1 "*'*; Twist2 Cre/+ mice (Fig. 5B). FGF 
signaling triggered the MAP kinase cascade, resulting in the 
phosphorylation and activation of p44/42 (Erkl/2) and the 
expression of two effector molecules, Erm and Pea3 [34]. 
Consistent with the reduced Fgf2 and Fgfr expression, the 
immunohistochemistry revealed that the levels of phospho-Erkl/ 
2 were considerably lower in TwistF 7 " x/+ ; Twist2 Cre/+ mice than in 
the control mice (Fig. 5C). Accordingly, quantitative real-time 
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Figure 6. Effects of Twistl, Twist2 and E12 on the activity of a 4.9 kb Fgfr2 promoter fragment. C3H10T1/2 (A) and MC3T3-E1 cells (B) 
were transiently co-transfected with a 4.9 kb Fgfr2 promoter luciferase construct and the indicated expression constructs, along with a pRL-TK 
construct as an internal control. The luciferase activities were determined by a dual luciferase assay system, and the promoter activities were 
expressed as luciferase activities relative to that of the control. The values represented mean ± SD. n = 3 for each group, "a" indicates significant 
difference from the control (p<0.05); "b" denotes a significant difference from all other groups (p<0.05). 
doi:10.1371/journal.pone.0099331.g006 



PCR demonstrated that the levels of Erm and Pea3 transcripts were 
significantly downregulated (Fig. 5A). These data demonstrated 
that FGF signaling was reduced in the TwistF 1 '"^; Twist2 Cre/+ 
mice, suggesting that Twistl and Twist2 might upregulate the 
expressions of Fgf2 and Fgfrs. 

Twistl and Twist2 stimulated Fgfr2 promoter activity 

As described above, Fgfr2 was significantly downregulated at 
both the mRNA and protein levels in the TwistP bx/+ ; Twist2 Cre/+ 
mice. The skeletal phenotype of TwistP lox/+ ; Twist2 Cre/+ mice 
partially resembles that of the Fgfr2 conditional knock-out mice 
[22]. Therefore, we examined whether Twistl and Twist2 could 
upregulate the Fgfr2 promoter activity in vitro. Our promoter 
luciferase reporter assay revealed that Twistl or Twist2 alone was 
unable to significantly stimulate a 4.9 kb Fgfr2 promoter fragment. 
However, they strongly enhanced the stimulatory activity of El 2, a 
ubiquitously expressed Twist binding partner, in both the 
C3H10T1/2 mesenchymal cells and MC3T3-E1 pre-osteoblast 
cells (Fig. 6). This in vitro evidence implied that Twistl and Twist2, 
together with El 2, could regulate the Fgfr2 expression. 

Discussion 

Twistl and Twist2 are two highly conserved members of the 
Twist subfamily of bHLH transcription factors. In this study, we 
generated Twistl- and TojirfJ?-haploinsufficient mice (TwistP'" x/+ ; 
Twist2 Cre/+ ), and found that these mice displayed delayed closure 
of the cranial sutures. Furthermore, we found that the Twist f° x/+ ; 
Twist2 Cre/+ mice presented with reduced bone formation and 
impaired osteoprogenitor proliferation and osteoblast differentia- 
tion as well as compromised FGF signaling. 

Previous studies have demonstrated that the cranial sutures of 
the Twistl heterozygous mice prematurely fuse [10,11,32]. To 
confirm this finding, we crossed the Twutl -&oxed mice {Twistl flox/ - /k ") 
with Sox2-Cre transgenic mice expressing a Cre recombinase 
ubiquitously in the epiblast cells at embryonic day 6.5 and 
generated conditional Twistl -haploinsufficient mice [TwistP lox/+ ; 
Sox2-Cre) (Figure S3). We found that the Twistf° x/+ ; Sox2-Cre 



mice developed craniosynostosis - a phenotype that is similar to 
the Twistl heterozygous mice [11,32], characterized by much 
narrower sagittal and interfrontal sutures and had additional toe 
originating from a duplicated metatarsal as demonstrated by 
alcian blue and alizarin red staining (Figure S3A-D). The plain 
X-radiography and histological examination did not reveal 
apparent difference in the radiopacity of the long bones, in the 
metaphyseal trabecular bones or in the diaphyseal cortical bones 
between TwistP 1 "" 1 ^; Sox2-Cre and control mice (Figure S3E-G). 
In contrast, although the Twist2 heterozygous mice are viable 
and fertile with no apparent phenotype [20,31], the TwistP hx/+ \ 
Twist2 C " /+ mice showed reduced growth of both cranial and 
long bones as well as delayed closure of the cranial sutures. 
Therefore, the TwistP 1 '"^; Twist2 Cre/+ mice displayed a skeletal 
phenotype that is opposite to that of Twistl heterozygous mice. 

Our current studies and the studies from other labs have 
demonstrated that Twistl and Twist2 share similar in vitro and in 
vivo functions [11]. Therefore, the phenotypic difference between 
single Twistl /J?-heterozygous mice and Twistl^" x/+ ; Twist2 Cre/+ 
mice might be due to the overlapping expression patterns of Twistl 
and Twist2 during mouse embryonic development [6] . Both Twistl 
and Twist2 are highly expressed in the condensed mesenchyme 
that later gives rise to osteoblasts and chondrocytes [6,11,22]. 
Thus, the outcomes of the loss of Twist 1 and/ or Twist2 really 
depend on the relative expression of both genes in a specific tissue. 

Accumulating evidence supports the notion that Twistl might 
control cranial suture development through modulating FGF 
signaling. It was found that the mutations in FGF receptors 
FGFR1, FGFR2, and FGFR3 in humans are associated with 
craniosynostosis, a characteristic phenotype of the Saethre- 
Chotzen Syndrome caused by dominant loss-of-function TWIST 1 
mutations [35]. In addition, the primary cranial osteoblasts 
isolated from SCS patients with Twistl mutations show reduced 
FGFR2 transcript levels, which can be restored by overexpression 
of TWIST1 [18]. Twistl haploinsufficiency in mice also results in 
an altered Fgfr2 expression pattern in the cranial sutures [17]. It 
has been proposed that Twistl haploinsufficiency favors the 
formation of Twist 1 homodimers in the osteogenic front of cranial 
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sutures, which results in the upregulation of Fgfr2 expression and 
leads to craniosynostosis [16]. Taken together, the data from these 
human and mouse studies suggest that Twist 1 might regulate FGF 
signaling, particularly the Fgfr2 expression, in a context-dependent 
manner. 

Our current studies provide further support that Twist 1 and 
Twist2 regulate FGF signaling in bone formation. We found that 
Twist f" x/+ ; Twtst2 c " /+ mice had reduced FGF signaling in bone as 
a result of decreased expression of not only Fgfr2 but also Fgfrl, 
Fgfr3, Fgfr4 and even Fg/2. FGF signaling plays essential roles 
throughout osteogenesis, including the commitment of mesenchy- 
mal cells to osteoprogenitors, the proliferation and differentiation 
of osteoprogenitors to osteoblasts, and osteoblast apoptosis [22,36]. 
Indeed, the Twistf" x/+ ; Twist2 Cre/+ mice displayed reduced 
proliferation of osteoprogenitor cells and defective osteoblast 
differentiation. Some of the skeletal abnormalities even resemble 
those of Fgfr2 conditional knockout mice [22]. Consistent with the 
in vivo results, our in vitro promoter luciferase assays further 
supported the role of Twist 1 and Twist2 in the upregulation of 
Fgfr2 promoter activity when E12 was present. However, it 
remains to be determined why only Twist/El 2 heterodimers, 
instead of Twist homodimers, stimulate Fgfr2 promoter activity. 

In summary, our current study suggested that Twist 1 and 
Twist2 may synergistically enhance the proliferation and differ- 
entiation of osteoprogenitors via the upregulation of FGF signaling 
during skeletal development. However, further studies are 
necessary to determine whether these two genes perform identical 
functions and could completely replace each other in vivo. 
Research is also needed to advance our understanding of how 
the expression levels of Twist 1 and Twist2 are temporally and 
spatially regulated during development and how fine-tuning 
achieves the optimal Twist protein level needed for normal 
skeletal development. Such knowledge is essential for the 
development of future therapies aimed at correcting the effects 
of Twist deficiency in humans. 

Supporting Information 

Figure SI The mRNA levels of Twist 1 and Tvuist2 in the 
long bones of TvuistP lox /+ ; Tvuist2 C ™ ' /+ mice. Real-time 
PGR was performed with total RNA isolated from the long bones 
of the 8-day-old control and TwistF 1 '"^; Twist2 c " /+ mice. The 
primers used for Twistl were sense 5'- CAGCGGGTCATGGC- 
TAAC-3' and antisense 5'- GCAGGACCTGGTACAGGAAG- 
3', and for Twist2 sense 5'- AGCAAGAAATCGAGCGAAGA-3' 
and antisense 5'- C AGCTTGAGC GTCTGGATCT-3 ' . The 
mRNA levels of Twistl and Twist2 were about three folds less in 



the Twistl flox/+ ; Twist2 Cre/+ mice than in the control mice. The data 

represented three analyses (« = 3) for each group. 

(TIF) 

Figure S2 Osteoblast apoptosis in the TvuistP 1 "" /+ ; 

Tmist2 Cre/+ mice. TUNEL assay was used to analyze the 
osteoblast apoptosis in the long bones of 6-day-old Twist^ bx/+ ; 
Twist2 Cre/+ mice and control mice. Three serial sections from each 
of four individual TwistF q " x/+ ; Twist2 c " /+ mice and control 
littermates were counted. No significant difference in osteoblast 
apoptosis was found between the two groups of mice. 
(TIF) 

Figure S3 Skeletal abnormalities of TvuistP l '" c/+ ; Sox2- 
Cre mice. (A) The skeletons of 7-day-old TwistV^" xl ^ (control) and 
TwistF 1 '"^; Sox2-Cre mice were stained with alcian blue and 
alizarin red. (B— C) Alcian blue and alizarin red stained skulls, 
femurs and tibiae, and hind feet from 7-day-old control mice and 
Twistf° x/+ ; Sox2-Cre mice; The Twistf" x/+ ; Sox2-Cre mice showed 
much narrower sagittal and interfrontal sutures (arrows; B) and 
had additional toe (arrow) originating from a duplicated metatarsal 
(MT; D), but the femurs and tibiae showed no apparent difference 
between two groups (C). (E) Plain X-radiography of the hind limbs 
from 7-day-old control and TwistF^" x, ^\ Sox2-Cre mice. No 
apparent difference was noted between the two groups of mice. 
(F and G) Histological examination of TwistP hx/+ ; Sox2-Cre mice. 
Tibia sections of 7-day-old control and TwistP 1 '"^; Sox2-Cre mice 
were stained with H&E. No apparent difference was observed in 
the metaphyseal trabecular bone (F) or in the diaphyseal cortical 
bone (G) between these two groups. 
(TIF) 

Materials and Methods SI Generation of TwistP 1 "" /+ ; 
Sox2-Cre mice. 

(DOCX) 
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